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Abstract. A 1-D loop radiative hydrodynamic model that incorporates the effects of gravitational stratification, 
heat conduction, radiative losses, external heat input, presence of helium, and Braginskii viscosity is used to 
simulate elementary flare loops. The physical parameters for the input are taken from observations of the Bastille- 
Day flare of 2000 July 14. The present analysis shows that: (a) The obtained maximum values of the electron 
density can be considerably higher (4.2 x 10^^ cm~^ or more) in the case of footpoint heating than in the case 
of apex heating (2.5 x 10^^ cm~^). (b) The average cooling time after the flare peak takes less time in the case 
of footpoint heating than in the case of apex heating, (c) The peak apex temperatures are significantly lower (by 
about 10 MK) for the case of footpoint heating than for apex heating (for the same average loop temperature 
of about 30 MK). This characteristic allows us to discriminate between different heating positioning, (d) In both 
cases (of apex and footpoint heating), the maximum obtained apex temperature T"^""^ is practically independent 
of the heating duration at, but scales directly with the heating rate Eho- (e) The maximum obtained densities 
at the loop apex, n™"^, increase with the heating rate Eho and heating duration at for both footpoint and apex 
heating. In Paper II we will use the outputs of these hydrodynamic simulations, which cover a wide range of the 
parameter space of heating rates and durations, as an input for forward-fitting of the multi-loop arcade of the 
Bastille-day flare. 
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parameters. Hydrodynamic modeling or nare loops, how- 

ever, have been performed for single flare loops (MariskaJ t,i- ■ r . .iiri- 
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tion studies have been orchestrated in a multi-loop con- 
figuration (Hori ct al. 1997, 1998). Even the multi-loop 
simulations have been designed only in the simplest way. 



fiares. In this Paper I we perform numerical simulations 
with a 1-D hydrodynamic code to obtain the temperature 
evolution r™°^(t), Tg'^''9(t) and density evolution n™"^(t) 
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tion (footpoint, apex) of the heating functions Enit). In 
Paper 11 we parameterize the results, suitable to forward- 
fitting of a multi-loop system and fit the multi-wavelength 
data of the Bastille-day flare 2000 July f 4 using TRACE, 
Yohkoh/SXT, HXT, and GOES data. In subsequent pa- 
pers we plan to extend the hydrodynamic results of this 
flare to constrain magnetic modeling, magnetic reconnec- 
tion geometries, and particle acceleration processes. 

2. The Radiative Hydrodynamic Model 

In order to describe plasma dynamics in a coronal loop, 
we solve numerically ID radiative hydrodynamic equa- 
tions (infinite magnetic field approximation) that resemble 
closely the Naval R esearch Labora tory Solar Flux Tube 
Model (SOLFTM) (< Mariskalll987i) . The numerical code 
that we use is a ID ver sion of the Lagrangian Re-map 
code l|Arber et al. Il200l[l with radiative loss limiters. As 
in SOLFTM the coronal loop is connected with a dense, 
cold (lOOOOK), 5 Mm thick plasma region, that mimics 
the chromosphere, and which because of its large density 
provides sufficient amount of matter to fill the loop during 
the flare. The use term of 'chromosphere' should be taken 
with caution, as we do not consider realistic one (with 
ionizations, proper radiative transfer, etc.). Instead, we 
simply mimic chromosphere by considering it as a dense, 
cold (IGOOOK) plasma, which provides source of matter 
inflow into the loop during the flare, while we use ra- 
diative loss limiters to prevent plasma from catastrophic 
cooling, as we use optically thin plasma approximation 
(ex rig). In such approach, specific initial conditions in re- 
gion connecting corona to the chromosphere and chromo- 
sphere itself have little or no effect on the corona dynamics 
as they are rapidly modified, self-consistently, according 
to the radiative hydrodynamics equations. In brief, our 
model includes: the effects of gravitational stratification, 
heat conduction, radiative losses, added external heat in- 
put, presence of helium, hydrodynamic non-linearity, and 
Braginskii bulk viscosity. For the radiative loss function 
we use the following parameterization. 



10-26.60 
10-17.73 

10-21.94 
10-10.40 rp-2 
10-21.2 
10-31 r^2 
10-21.85 
10-48.31 rpe.l5 
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T > 10^-6 

106.3 < lO^-s 

105.8 < y < 106.3 

105.4 ^ 105.8 

lO^-^ < T < lO^-* 
104.6 < lO^-^ 
lO^-^ < T < lO^-^ 

103.9 <:T < lO^-^ 

103.6 < 10^-^ 



which is an extens ion oflRos ner et al. | l|l978 [) com piled 
from other sources l|Priest lll982llPeres et al. Ill982() . 



3. Numerical results 

We start numerical simulations from the following configu- 
ration: we take a semicircular loop with a length of L = 55 
Mm (which corresponds to an average loop arcade radius 



of r = 17.5 Mm, as derived bv lAschwanden fc Alexander I 
l|200l|) . Table III, for the Bastille-day flare). We keep the 
coronal part of the loop initially at a temperature of 1 MK 
and at a mass density of p = fimpUe = 6.6 x 10~^^ g cm~^ 
(at the loop apex), for a helium-to- hydrogen number den- 
sity ratio of 0.05, i.e., with a mean molecular weight of 
fj.= {1 + 0.05 X 4)/(l + 0.05 X 2) = 1.1, this corresponds 
to an electron density of Ue ~ 3.6 x 10® cm~'^. The res- 
olution in all our numerical runs was fixed to 1000 grid 
points, which were distributed non- uniformly in order to 
properly resolve strong gradients in the region connecting 
corona to the chromosphere. As a convergence test, runs 
with 3000 grid points were made, which showed no dif- 
ference to the case with 1000 grid points, thus confirming 
sufficient numerical resolution in our simulations. 

3.1. Heating Function 

The heating function in flare loops is probably quite dif- 
ferent from that of non-flaring loops. Flare loops are 
filled by up-flowing heated plasma from the chromo- 
sphere, once the chromospheric footpoints become im- 
pulsively heated from precipitating non-thermal particles 
and/or downward propagating hot thermal conduction 
fronts from the coronal reconnection site. Chromospheric 
evaporation seems to be the main matter inflow source 
for flare loops. The heating function has therefore to 
accommodate very localized heati ng at the footpoint s 
l|Tsiklauri fc Nakariakov ■■200L.Aschwanden et al. Il200l|) . 
On the other side, reconnection outflows from the recon- 
nection region contain also heated plasma and heat the 
flare loops from the ap ex side, for instance in the s tandard 
reconnection model of IKodd fc Pneuman I l|l976l) . which 
seems to fit the magn etic configuration of the Bastille - 
day flare to first order jAschwanden fc Alexander1l200lf) . 
Thus, to allow for both options, we have performed numer- 
ical simulations for both cases separately, i.e., for apex and 
footpoint heating functions. 

We have used the following heating function in our 
simulations: 
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Here, Efj{s) and E\j{t) are the spatial and temporal 
parts of the heating function, which are taken to be in- 
dependent of each other for simplicity. Eq is the heating 
rate in units of [erg cm~'^ s~^]. The positions s = ±so 
are the locations with the maximum heat deposition, i.e.. 
So = for apex heating. The heat deposition length scale 
is called (i.e., the spatial width of the Gaussian). The 
te mporal part of the heating function is similar to one used 
bv lAschwanden fc Alexander I l)200l[l (cf. their Eq. (31)), 
where tp is the flare peak time and at is the duration of the 
flare (i.e., the temporal width of the Gaussian). However, 
our choice of the temporal part of the heating function 
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is such that there is a small background heating present 
at all times (either at footpoints or the loop apex) which 
ensures that in the absence of flare heating (when a = 0, 
a parameter that determines the flare heating amplitude) , 
the average loop temperature stays at 1 MK. For easy 
comparison between apex and footpoint heating cases we 
fix the flare heating amplitude Qp at a given different value 
in each case. This ensures that the average loop tempera- 
ture peaks at about at the observed value of 30 MK in both 
cases, when the the flare heating is on (a = 1). Then we 
vary also a = 0.01,0.05,0.25,0.5,1.0 to obtain different 
apex temperatures T™"^. Hereafter, we denote the (flare) 
heating rate by Ehq — EoaQp. 

In all our numerical runs presented here, 1/(2(Ts) was 
fixed to a value of 0.01 Mm~^, which gives a heat depo- 
sition length scale of as — 7 Mm. This is a typical valu e 
determined from observations ijAschwanden et al. lEool . 
The flare peak time was fixed in our numerical simulations 
at tp — 822 s. The time step of data visualization (which 
in fact is much larger that the actual time step, 0.035 s, 
in the numerical code) was chosen to be At = 10 s. 

3.2. Case of Apex Heating 

In the case of apex heating we fix sq = Mm in Eq.(l). 
Initially we run our numerical code without flare heat- 
ing, i.e. we put a = (in this manner we turn off flare 
heating). In all our numerical runs we keep background 
heating on at all times, so that flare occurs starting from 
a steady, equilibrium loop. The result of this simulation is 
presented in Fig. 1. Dash-dotted curves in panels (a) and 
(b) show the evolution of the average and apex tempera- 
tures in time. Since the observational data of the Bastille 
day flare, namely the temperature evolution (cf. to p left 
panel of Fig. 11 in lAschwanden fc Alexander I l)200l|) '). was 
derived for the entire fleld of view of the instrument, we 
found that it is useful, in addition to the apex tempera- 
ture, to track also the average loop temperature, defined 
as the sum of temperatures along the loop divided by the 
number of grid points at a given snapshot. Since the out- 
puts from our hydrodynamic simulations will be used to 
model the observational data of the Bastille-day flare in 
paper II, Eq was fixed to 0.002 erg cni^"^ s^^ in all runs of 
this subsection, which insures that the average loop tem- 
perature stays at 1 MK in the absence of flare heating 
(cf. Fig. la). Note that the apex temperature in this case 
tends to have a higher asymptotic value of 2.5 MK (cf. 
Fig. lb). The obtained asymptotic value of the loop apex 
number density is = 0.9 x 10^ cm^'^ in this case (and 
thus is not visible in Fig. Ic). 

Then, we run our numerical code with flare heating, 
i.e. we put a ~ 1 and fix Qp = 3 x 10*, so that it yields 
peak average temperature of about 30 MK, as observed 
during the Bastille-day flare ijAschwanden k, Alexander! 
L2OOI), and then we vary the duration of the heating phase 
to at = 41, 164, 329 s (i.e., the temporal width of the 
Gaussian). Note that it is the combined action of the 



heating rate Qp, the flare heating amplitude a, and the 
duration at of the flare, that determine the effect of flare 
heating, i.e. the term aQpCxp (— (t — tp)^ /2af) in Eq.(l), 
compared with the background heating rate (a = 0). The 
results are presented in Fig. 1. There are several notewor- 
thy features in this graph: (1) As expected, an increase of 
the heating time interval at yields an increase of the flare 
duration (cf. Fig. l(a,b)). (2) However, an increase in the 
duration of the heating phase does not affect the obtained 
maximum values of the average nor the apex temperature 
(cf. Fig. la and lb). This is counter-intuitive, because we 
expect that the amount of deposited heat increases with 
the heating duration. This invariance of the obtained flare 
temperatures with respect to duration of the flare could 
probably be explained by some form of saturation (or bal- 
ance) in the combined action of losses in the system - heat 
conduction and radiative losses - and flare heat input. (3) 
An increase in the duration of heating should naturally re- 
sult in an increase of the plasma density in the loop. For 
the long duration flares (the thickest solid line in Fig. Ic) 
there are three clearly different physical regimes which yet 
have to be properly identified. (4) Very useful information 
can be extracted from the velocity outputs at the foot- 
points (s = ±27.5) and at the apex (s = 0), as function 
of time (cf. Fig. Id, le, If). We find that strong up-flows 
(up to 400 km s"""^) are present at the footpoints of the 
flare loops during the flare onset, as the deposited heat 
(delivered by conduction from the apex to the footpoints) 
causes material evaporation from the dense chromosphere. 
The late up-flow phase is followed by an oscillatory phase 
with typical amplitudes of a few tens of km s^^, which in 
turn is followed by down-flows (up to 100 km s^^), when 
plasma is drained out of the loop (5) Note that the ve- 
locity dynamics at opposite footpoints remains perfectly 
symmetric at all times (cf. Fig. Id, le. If), while net flow 
through the apex remains zero at all times as expected. 
This is due to the symmetry of the problem. (6) The peak 
up-flow velocity during the onset of a flare decreases with 
the increase of the duration of the flare (for at = 41s). 
In Fig. Id the peak up-flow velocity is about 400 km s^^, 
while for at — 164s (Fig. le), it is about 200 km s^^). 
This may seem counter-intuitive at first glance, since we 
would expect that the duration of the up-flow itself in- 
creases with increasing heating duration at. Thus, the net 
material evaporated into the loop from chromosphere still 
increases with the increase of the heating duration, so that 
the densities obtained at the apex increase (cf. Fig. Ic). 
(7) Yet another interesting observation is that the time in- 
stances when down-flows abruptly end after the flare (cf. 
Fig. Id, le. If), these times correspond exactly to the same 
time instances when both the average and apex temper- 
ature curves reach these asymptotic values (thin, thick, 
and very thick solid lines that join the dash-dotted lines 
in Figs, la and lb), signaling the end of the flare heating 
phase and the onset of a steady state. In fact, statements 
(6) and (7) are true for all numerical runs performed (see 
below) . 
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Fig. 1. (a) Average temperature as a function of time in the case of apex heating. Dash-dotted curve corresponds to 

the case of background heating (a = 0). Tlie tliin. tliick and thickest solid lines correspond to the cases of background 
plus flare heating (a = 1) for different flare durations, at = 41, 164, 329 s, respectively, (b) The same as in (a), but for 
the temperature at the apex, (c) The same as in (a), but for the number density at the apex, (d) Velocities versus 
time at footpoints (s = ±27.5), solid and dotted curves, and apex (s = 0), dashed curve, for the case of background 
plus flare heating (a = 1) for the flare duration, CTj = 41 s. (e) The same as in (d), but for = 164 s. (f) The same 
as in (d) , but for at = 329 s. 



In a next step we investigate the effect of the flare 
heating amplitude, by fixing the duration of the flare 
fTt = 41 s and varying a = 0.01, 0.25, 1.0 (again with fixed 
Qp = 3 X 10^). The results of these numerical runs are 
presented in Fig. 2. There are several interesting features 
in this graph: (1) As expected, a decrease of a results in 
a decrease of the obtained flare temperature (cf. Fig. 2a 
and 2b). (2) Also, a decrease of the flare heating ampli- 
tude does not affect the duration of the flare (cf. Fig. 2a 
and 2b). (3) In addition, a decrease of the flare heating 
amplitude naturally results in a decrease of the plasma 
density in the loop. This is understandable as the less 
deposited heating rate causes less material evaporation 
form the chromosphere and, hence, produces less dense 
and cooler loops during the flare (cf. Fig. 2c, 2a, and 2b). 
(4) Again, as in Fig. Id, le. If, we studied the resulting 
velocities at the footpoints (s = ±27.5 Mm) and apex 
{s = 0) as function of time (cf. Fig. 2d, 2e, 2f). Our sim- 
ulations show that the peak up-flow velocities during the 
onset of the flare decrease with decreasing heating ampli- 



tudes. For Q = 1.0 (Fig. Id), the peak up-flow velocity 
is about 400 km s~^, while for a = 0.01 (Fig. If), it is 
about 160 km s~^). Also, note that the 4-fold decrease in 
the heating flare amplitude (Fig. le) still results in down- 
flows of the order of 100 km s~^, which are typical for 
most runs, while further decrease in a (Fig. Id) causes an 
absence of any noticeable down-flows. 

3.3. Case of Footpoint Heating 

In the case of footpoint heating we fix sq = ±30 Mm in 
Eq.(l), i.e. the (spatial) peaks of the heating are chosen to 
be at the bottom of the region connecting corona to the 
chromosphere (i.e. top of chromosphere). As for the case 
of apex heating, we run initially our numeric code without 
flare heating (a = 0). The output from this numerical run 
of the code is presented in Fig. 3. Eq is in all runs in this 
subsection fixed to a value of 0.01 erg cm~^ s~^, which 
insures that the average loop temperature stays at 1 MK 
in the absence of flare heating (cf. Fig. 3a). Note that in 



Tsiklauri et al.: Radiative HD Modeling of solar flares 



5 



Average Temperature 



Apex Temperature 



Loop Apex N. Density 




1000 2000 3000 4000 
Time (s) 

Velocity at Footpoints and Apex 





1000 2000 3000 4000 
Time (s) 

Velocity at Footpoints and Apex 

400 

200 





1000 2000 3000 4000 
Time (s) 

Velocity at Footpoints and Apex 

400 

200 - 



(f) 




1000 2000 3000 
Time (a) 



1000 2000 3000 
Time (s) 



1000 2000 3000 
Time (s) 



Fig. 2. (a) Average temperature as a function of time in the case of apex heating. Dash-dotted curve corresponds to 
the case of background heating (a = 0). The thin, tliick and thickest solid hnes correspond to the cases of background 
plus flare heating (a = 1) for different flare peak amphtudes, a = 0.01,0.25,1.0, respectively (these correspond to 
flare heating rates Ehq = 0.6, 15.0, 60.0 erg cm~^ (b) The same as in (a), but for the temperature at the apex, 
(c) The same as in (a), but for the number density at the apex, (d) Velocities versus time at footpoints (s = ±27.5), 
solid and dotted curves, and apex (s = 0), dashed curve, for the case of background plus flare heating (flare duration 
at = 41 s) for the flare peak amplitude, a = 1.0 (which corresponds to flare heating rate Ehq = 60.0 erg cm~^ s~^). 
(e) The same as in (d), but for a = 0.25 {Ehq = 15.0 erg cm~^ s~^). (f) The same as in (d), but for a = 0.01 
{Eho = 0.6 erg cm~^ s~^). 



the case footpoint heating, the amount of deposited heat 
required to keep the loop at an average temperature of 1 
MK is 5 times greater than in the case of apex heating 
(recall that in the later case Eq = 0.002 erg cm~^ 
This increased heating is required in order for the upward 
heat conduction to be sufficient enough to keep the loop 
at this temperature. The apex temperature in the case 
of footpoint heating approaches an asymptotic value of 
1.7 MK (cf. Fig. 3a), while the obtained asymptotic value 
of the loop apex number density is 0.02 x 10^^ cm~^ in 
this case (not visible in Fig. 3c). Note that the asymptotic 
value of the loop apex temperature, T™"^ = 1.7 MK, is 
lower in the case of footpoint heating (due to the flatness of 
the temperature spatial profile along the loop) compared 
with the case of apex heating (T™"^ = 2.5 MK). The 
larger amount of deposited heat in the case of footpoint 
heating results in an asymptotic value of the the loop apex 
density (rie = 0.02 x 10" cm'^) that is more than twice 



the value for apex heating (rig = 0.009 x 10^^ cm~^). This 
is due to the fact that the more depositcxi hc^at causes more 
material evaporation from the chromosphere into the loop. 

Yet another interesting point that only appears in the 
case of footpoint heating is that there is a minimal heat 
deposition length scale, Us, which allows the existence of 
a steady loop at an average temperature of 1 MK. The 
results presented in this subsection are for a-g = 7 Mm. 
Smaller values cause a gradual formation of a condensa- 
tion at the apex of the loop (clearly seen in the animations, 
not included here) , which then forms a prominence due to 
the thermal instability and causes an eventual disappear- 
ance of the loop. This occurs even with a substantial in- 
crease of the heating amplitude Eq . This is caused by the 
inefflciency of upward heat conduction to keep the loop at 
a typical coronal temperature of the order of T « 1 MK. 

Our next step is to switch on flare heating as described 
in previous subsection. The physical parameters used here 
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Fig. 3. The same as in Fig. 1, but for footpoint heating. 



are also the same (apart from differences in the steady 
state as described above) for easy comparison. Namely, 
we put a = 1, and fix Qp = 1.5 x 10*, so that it yields an 
average peak temperat ure of about T™"^ = 30 MK as i n 
the Bastille-day flare ijAschwanden fc Alexander! l200l|) . 
and we vary duration of the flare by changing the heat- 
ing length scale to at = 41, 164, 329 s. The results of this 
simulation are presented in Fig. 3. Since, the dynamics of 
flare with variation of its duration has been described in 
considerable detail for the case of apex heating (Fig. 1), 
we focus now our attention on differences between the two 
cases. These can be summarized as follows: (1) We observe 
(Fig. 3c) that the obtained maximum values of the den- 
sity are considerably higher (ne = 4.2 x 10'^"'^ cm~'^) in the 
case of footpoint heating than in the case of apex heating 
(ne = 2.5x10^^ cm~'^), (compare Fig. 3c with Fig. Ic). This 
is due to the fact that footpoint heating is more efficient in 
evaporating material from the region connecting corona to 
the chromosphere and chromosphere itself, yielding denser 
loops during the flare. In fact, previous hydrodynamic sim- 
ulations assert that there is a maximum density limit, no 
matter how much heat is deposited. This limiti ng density 
is about 2 x 10" c m"^ ( llMariska et al. Ill989t) . cf. their 
Fig. 2, after 60 sV iFisher & HawlevI l(l99n>) get similar 
values. An analogo u s pro blem exists in quite Sun loops. 
lAschwanden et ahl l)200ll) points out that the observed 



TRACE loops show a higher density and pressure than ex- 
pected from the RTV law, which we might call the "over- 
density loop problem" . However, as we can see in Fig. 3c, it 
is possible to break that "density limit" and the reason for 
this is that we have used footpoint heating. In the case of 
footpoint heating it is possible to evaporate more material 
from the chromosphere and the region connecting corona 
to the chromosphere, where the density is orders of magni- 
tude higher than in the corona. In the case of apex heat- 
ing, in contrast, insufficient downward h eat conduction 
prevents significant evaporation. Works of iMariska et al. 1 
(1989) and Fisher & Hawlcv ( 1990) presumably were mo- 
tivated to use a pex heating because of the standard recon- 
nection model |)Kopp fc Pneuman Ill976|) . which implies 
that t he heating funct i on is l ocalized in the loop apex (ac- 
tually, iMariska et al. I l)l989|) used a model in which the 
(coUisional) heating by a beam of nonthermal electrons 
injected at the apex followed the density evolution of the 
flare loop and so varied in time in response to the evolu- 
tion of the flare). Consequently they obtain lower densities 
in the loop. Similar results that by modeling a loop as set 
of small scale, impulsively heated filaments one can gen- 
er a.llv reproduce t he observed flaring loops were obtained 
bv lWarrenet all j2003). (2) In the case of footpoint heat- 
ing, a flare takes less time on average (compare Fig. 3a 
and 3b with Fig. la and lb) - so the cooling of the loop 
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happens faster. This can be explained by the following 
reason: At the initial stages of cooling conductive losses 
dominate over radiative losses (when T > 20 MK), while 
the situation reverses (when T < 20MK) as the loop cools 
down and becomes less dense - plasma is flowing out of 
the loop (strong down flows were seen to accompany the 
cooling process in all of our numerical simulations and 
this is also confirmed in a number of observations). The 
time scale of conduction loss is proportional to the den- 
sity, while the time sca le of radiative loss is pro portional 
to the inverse density IjAschwanden fc Alexander .20011. 
Therefore, since the radiative losses dominate over the 
heat conduction losses for most of the time, it is natural 
to conclude that the denser loops would cool faster. (3) In 
the case of footpoint heating the peak apex temperatures 
{corresponding to the same, as in the case of apex heat- 
ing, average temperature of about 30 MK) are significantly 
lower (less by about 10 MK). The observational data of 
the Bastille-day flare, namely the temperature evolution 
(cf. to p left panel of Fig. 11 in lAschwanden fc Alexander I 
l|200lh V was derived for the entire field of view of the 
instrument, i.e. it tracks the dynamics of the average tem- 
perature. However, if one would have additional tempera- 
ture dynamics at a loop given point such as at the apex, 
our simulations would allow to discriminate between the 
different heating functions of the loop during the flare. 
(4) In the case of footpoint heating up-flow velocities are 
somewhat higher (compare Fig. 3d, 3e, 3f with Fig. Id, le. 
If). 



To complete the comparison between apex and foot- 
point heating cases, we investigate the effect of the flare 
heating amplitude, by fixing the duration of the flare to 
(Jt = 41 s, and vary a = 0.01,0.25, 1.0 (again with fixed 
Qp = 1.5 X 10''). The results of these numerical runs are 
presented in Fig. 4 (see for comparison Fig. 2). Most of 
the features seen in the case of apex heating are still valid 
in this case. A notable difference is that, in the considered 
parameter space, spikes in the apex density time profiles 
occur (Fig. 4c), which can be also seen in Fig. 3c after 
the flare. In fact, these spikes are present in all runs with 
footpoint heating. They occur when dense blobs of plasma 
(prominences), formed by the thermal instability, swipe 
through the apex. One should mention, however, that a 
ID code does not provide a fully adequate description of 
prominences, since it ignores the finite magnetic field ten- 
sion (we use an infinite magnetic field limit in our code). 
In fact, this is the reason why they do not stay steady at 
the apex, as the magnetic field in our case cannot bend to 
provide cavity for a stationary prominence. Yet another 
difference between the two cases are the higher obtained 
densities and up-fiow velocities during the flare in the case 
footpoint heating (due to more efficient evaporation) and 
lower obtained apex temperatures (due to the spatial flat- 
ness of the heating function along the loop) . 



3.4. Parametric Study 

In this subsection we present a parametric study of the 
problem by investigating the maximum obtained temper- 
atures and densities at the loop apex as a function of 
flare duration, at and the heating rate Eho — E^aQp. We 
have performed simulations for at = 41, 82, 164, 329 s and 
a = 0.01, 0.05, 0.25, 0.5, 1.0 with Qp = 3 x 10* for the case 
of apex heating and Qp = 1.5 x 10* in the case of footpoint 
heating. With appropriate multiplication by Eq in each 
case (i?o = 0.002 erg cm^^ s^^ for apex heating and Eq = 
0.01 erg cm~^ s~^ for footpoint heating), the fiare heating 
rates, Eho, are, 0.6,3.0,15.0,30.0,60.0 erg cm'^ s'^ in 
the case of apex heating and 1.5,7.5,37.5,75.0,150.0 erg 
cm~'^ s~^ in the case of footpoint heating. The results of 
these numerical runs are presented in Fig. 5. The actual 
data that was used to produce Fig. 5 is given in Table 1. 
In Figs. 5a and 5c we plot the maximum obtained tem- 
peratures, T™"^ (MK), at the loop apex in the case of 
apex and footpoint heating, respectively. We gather from 
these plots that: (1) In both cases T™"^ is practically in- 
dependent of the fiare duration at and it increases with 
the increase of the flare heating rate Ehq. (2) T™°^ is 
up to 10 MK higher in the case of apex heating than in 
the case of footpoint heating. This due to the steepness of 
the temperature profile along the loop in the case of apex 
heating (corresponding to the same average temperature 
as in the case of footpoint heating). In Figs. 5b and 5d we 
plot the maximum obtained densities, n^^^/lO^^ (cm^'^), 
at the loop apex in the case of apex and footpoint heating, 
respectively. It can be seen that (1) In both cases n™**^ in- 
creases with the increase of both the flare heating rate Ehq 
and flare duration at, and (2) n™"^ is about twice as high 
in the case of footpoint heating than in the case of apex 
heating. This is due to more efficient material evaporation 
from the chromosphere and the region connecting corona 
to the chromosphere in the case of footpoint heating. We 
have also performed parametric study varying the loop 
length. The results are presented in Table 2. Note that for 
consistent comparison with other numerical runs, in the 
case of footpoint heating, when varying loop length, we 
have shifted sq so that spatial maximum of the heating 
function always stayed at the top of chromosphere. 



4. Conclusions 

In summary, we have used a radiative hydrodynamic nu- 
merical code to simulate flares. The physical parameters 
of the input were obtained from observations of the the 
Bastille-day flare (Aschwandcn fc Al exander 200ll). Our 
simulations confirm the general picture of flare dynamics: 
Transient heat deposition either at loop footpoints (with 
maximum heat input at the bottom of the region con- 
necting corona to the chromosphere or i.e. top of chromo- 
sphere) or at the apex leads to an average loop tempera- 
ture of Tavg ~ 30 MK first. Then, evaporation of material 
from the chromosphere and the region connecting corona 
to the chromosphere into corona ensues with up-flows in 
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Fig. 4. The same as in Fig. 2, but for footpoint heating. However, the flare peak amplitudes a = 0.01,0.25,1.0 
correspond here to the flare heating rates Ehq = 1.5,37.5, 150.0 erg cm~^ ^^'^)- 



the order of a few hundreds of km s~^. During the peak of 
the flare, the combined action of heat input and conduc- 
tive and radiate loss yields an oscillatory flow pattern with 
typical amplitudes of up to few tens of km s~^. Finally we 
enter a cooling phase, when down-flows in the order of 
hundred km s^^ can be seen as the plasma drains out of 
the loop, ultimately reaching an equilibrium. 
We have established the following: 

1. In the case of footpoint heating, the obtained max- 
imum values of the density are considerably higher 
(4.2 X 10^^ cm~'^ or more) than in the case of apex 
heating (2.5 x 10^-^ cm"'^). This is due to the fact that 
footpoint heating is more efficient in evaporating ma- 
terial from the region connecting corona to the chro- 
mosphere and chromosphere itself, which yields denser 
loops during the flare. In the case of apex heating, 
which was used to model flares, insufficient downward 
heating conduction prevents significant material evap- 
oration. 

2. In the case of footpoint heating, as compared to the 
apex heating, on average cooling after the flare takes 
less time. This due to the fact that the time scale of 
conduction loss is proportional to the density, while the 
time scale of radiative loss is reciprocal to the density 



(lAschwanden &: Alexander Therefore, since ra- 

diative losses dominate over heat conduction losses for 
most of the time, it is clear that the denser loops cool 
faster. 

3. In principle, our simulations would allow to discrim- 
inate between different heating functions of the loop 
during the flare, if one would have temperature dy- 
namics in a given point of the loop, such as at the 
apex. This is based on our observation that in the case 
of footpoint heating the peak apex temperatures (cor- 
responding to the same, as in the case of apex heating, 
average temperature of about 30 MK) are signiflcantly 
lower (less by about 10 MK). 

4. In the case of footpoint heating, up-flow velocities are 
higher (roughly up to 100 km s""'^) than in the case of 
apex heating due to more efflcient evaporation. 

5. In both cases (of apex and footpoint heating) the max- 
imum obtained temperature T™"^ at the loop apex is 
practically independent of the heating duration at, but 
it increases with higher heating rates Eho- 

6. The maximum obtained densities at the loop apex in- 
crease with the increase of both the flare heating rate 
Eho and the heating duration at, in the case of apex 
as well as footpoint heating. 
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Fig. 5. (a) The dependence of maximum obtained temperatures at the loop apex as a function of flare duration, at 
and heating rate Eho for the case of apex heating, (b) The same as in (a), but for the maximum obtained density in 
the case of apex heating, (c) The same as in (a), but for the case of footpoint heating, (d) the same as in (b), but 
for footpoint heating. 



7. Varying the loop length (see Table 2) in the range 
of L = (0.25, ...,2.0) X Lo (with Lq = 55 Mm), we 
find (1) that the mean loop temperature averaged over 
the loop length does not change dramatically, (2) that 
the loop apex temperature increases notably for longer 
loops only for the case of apex heating, but much less 
for footpoint heating, and (3) that the mean electron 
density decreases somewhat with longer loops, i.e., 
rie/lO" (cm-3) = 10.45L^°-362 for apex heating and 
rie/lO" (cm-3) = 9.62L-°-207 foj. footpoint heating. 
Here, L is loop length in Mm. 

In practically all of our numerical runs we have de- 
tected quasi-periodic oscillations in all physical quanti- 
ties. In fact, suc h oscillations are frequ ently seen during 
solar flares (e. g. iTerekhov et al. I ll2002l)') as well as stel- 
lar flares fe.g. iMathioudakis et al. I l)2003j) '). Our prelimi- 
nary analysis shows that quasi-periodic oscillations seen 
in our numerical simulations bear many similar features 
as the observed ones. The key point is that the tradi- 
tional explanation of these oscillations in the observations 
involves MHD waves. In the numerical simulations pre- 
sented here, however, they are likely to be produced by 



standing sound waves caused by impulsive and localized 
heating. Therefore, our explanation of these oscillations 
is purely hydrodynamic - they are related to the stand- 
ing slow r node acoustic wave s, similar to the observed by 
SUMER JWang et al. II2002I) . A detailed study of these 
quasi-periodic oscillations will be presented elsewhere. 

In a next step we plan to use the outputs of this para- 
metric study of hydrodynamic simulations, which cover a 
wide parameter range of heating rates and heating time 
scales, as input for forward-fitting to the observed physical 
parameters (densities and temperatures) of the multi-loop 
flare on Bastille-day 2000 July 14. 
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Table 1. Flare loop temperatures and densities as function of the heating rate and heating duration, for a fixed loop 
length of L = 55 Mm, obtained from ID hydrodynamic simulations. Units of Eho are in (erg cm~^ s^^). at is measured 
in seconds, while temperatures and densities are given in (MK) and (cm~^) respectively. 
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Table 2. Flare loop temperatures and densities as function of the loop length (in Mm), for a fixed heating duration 
of at = 329 s. Units of Eho are in (erg cm"'^ s~^). Ut is measured in seconds, while temperatures and densities are 
given in (MK) and (cm~^) respectively. 
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